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CLASSES OF SMOOTH SOLUTIONS TO THE 
MULTIDIMENSIONAL BALANCE LAWS OF GAS DYNAMIC 
TYPE ON THE RIEMANNIAN MANIFOLDS 

OLGA S. ROZANOVA 


Abstract. The paper is devoted to the special classes of solutions to 
multidimensional balance laws of gas dynamic type. In the velocity field 
for the solutions of such class the time and space variables are separated. 
The simplest case is the solution with the linear proHle of velocity in the 
Euclidean space. 


1. Balance laws on the riemannian manifold 


1.1. Preliminaries. Let S be a smooth n - dimensional riemannian mani¬ 
fold with local curvilinear orthogonal coordinates x*, i = Consider the 

following system of nonlinear partial differential equations 

(ff(U,x,t)); + V,T 7 '^'(U,x,t) + Q*(U,x,t) =0, i,j = l,...,n. (1.1) 

Here U = (tt^,..., u"^), R = {R^, ..., R^) and Q = {Q^, ..., Q^) are the 
vector-functions on the tangent bundle of S [u is unknown), is a smooth 
tensor field, is a divergency of tensor. All differential operations with 

respect to space variables are performed in this coordinate system by means 
of smooth metric tensor gij of the manifold S, the variable t denotes the time, 
as usual. Recall that the covariant derivative of any contravariant vector 
X{X^,...,X^) is the tensor VjX^ = g + i,j,k = l,...,n, where 

are the Christoffel symbols. 

If the tensor field is sufficiently smooth, then the system (1.1) can be 
re-written in the form 


G;.(U, X, t)iu% + Ff (U, X, t)VkU^ + Q*(U, X, t) = 0, i = 
with new tensor fields G*- and or 


A0(U, X, t)^ + A^(U, X, t)|^ + Q(U, X, t) = 0, 

k=l ® 


n. 


( 1 . 2 ) 


(1.3) 


where = G), = Fj^, Q* = Q* + 

Remark that the system (1.2) includes only covariant derivatives, which 
coincide with the usual partial ones only in the case of the euclidean space. 


Date: July 11, 2002 . 


1 




2 


OLGA S. ROZANOVA 


Nevertheless, the representation of the covariant derivatives through the par¬ 
tial ones and the Christoffel symbols [ 2 ] adds the derivatives free terms only. 

We restrict ourself by the systems of gas dynamic type, namely, by the 
system of the Euler equation with the right hand sides of special form. 

The important fact is the possibility to reduce these balance laws to the 
symmetric hyperbolic form. 

Recall that the system of form (1.3) for the unknown function U(t,x) 
is called symmetric hyperbolic if the matrices ^■^(U,x,t),j = 0 , ...,n, are 
symmetric and the matrix is positively defined in addition. The 

following classical result is well known for the euclidean space M- ([3], i, 
particular cases in 0, 0, 0, 0): if the matrices A^{p,x,t) and the right 
hand side Q(p, x, t) smoothly depend on their arguments, have the continuous 
bounded derivatives up to order m +1 with respect to the variables (p, x) under 
bounded p, and the initial data U(0,x) and the function Q(0,x, t) belong to 
the Sobolev class m > l+n/ 2 , then the corresponding Cauchy prob¬ 

lem has locally in time unique solution in the class n^gC-^([0, T); 

Moreover, 

lim sup(||U|| r oo + llVxUll J.OO ) = -|-oo. 
t^T-0 

(Remark that on the initial data less strict requirements may be imposed, for 
example, the requirement of belonging uniformly local in space to 

[H], 0). 

As usual in the problem arising from physics the coefficients of system (1.3) 
depend only on the solution U. The fact simplifies signihcantly the formulation 
of the result. Namely, for the local in time existence of the Cauchy problem 
in the described class it is sufficiently to require the implementation of the 
condition Q(0) = 0 besides of the smoothness of coefficients. 

If the manifold S is diffeomorfic to 7^", then the theorem on the existence 
of smooth local in time solution can be transmitted to the case. 

The author don’t know the analogous result for the case of an arbitrary 
riemannian manifold. 

1.2. The system of the Euler equations. Eurther we shall consider the 
following system given on the riemannian manifold S in n dimensions with 
the metric tensor gij (below we denote it (E1-E3)) 

p(5tV + (V,V)V) = -Vp + pF(p,p,V,x,t), (1.2.1)(E1) 

dtp + div {p\) = Q ( 1 . 2 . 2 )(E 2 ) 

dtp + {y ,yp) + ^pdiv\ = Q. (1.2.3) (E3) 

Here p, p, V are the density, the pressure and the velocity vector, respectively. 

The constant 7 is the heat ratio, 7 > 1, moreover, in the physical case 7 < 

1 + ^. 
n 

The right-hand side pF is the forcing term, it is supposed to be smooth 
with respect to all arguments. 
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For example, it may describe the friction if F = with the non¬ 

negative function /i(f,x) and the constant a > 0, or the Coriolis term. In the 
last case F = [V x w], a; is a constant vector, corresponding to the vector of 
angular velocity (for n = 3) or F(V) = /V_i_, I = /(x, t) is the Coriolis parame¬ 
ter, VI = CijV^, Cij is the skew-symmetric discriminant tensor of the surfaces 
(for n = 2). 


2. The Euclidean space. Solutions with the linear profile of 

VELOCITY AND THEIR PROPERTIES 

We consider the classical solutions to (E1-E3) with the density sufficiently 
quickly decreasing as |x| —oo to guarantee a convergency of the integral 
f pixpdx (so called solutions with a finite moment). 

As well known, solutions to the Cauchy problem for system (E1-E3) may 
lose the initial smoothness for a finite time. Sometimes there is a possibility 
to estimate the time of a singularity formation from above (see, f.e., |10] and 
references therein). Moreover, in the case F = 0 the singularity appears in 
any solution with compactly supported initial data (f.e., [TT]b 

At the same time, there exist some nontrivial classes of globally smooth 
solutions. In the section our special interest will be the solutions with linear 
profile of velocity 

V = A(t)r + b(t), (2.1) 

where A{t) and b(t) are a matrix (n x n) and an n - vector, dependent on 
time, r is a radius-vector of the point. 


2.1. Construction of solutions with linear profile of velocity and in¬ 
tegral functionals. 


2.1.1. The system on the plane (n = 2). Now we construct solutions with 
linear profile of velocity for the important case F = LV, with matrix L = 

, /i = const >0,1 = const. This right-hand side describes in a 


/ -I 

V I -T 


simplest way the Coriolis force and the Rayleigh friction in the meteorological 
model neglecting the vertical processes. 

We consider only smooth solutions to (E1~E3) with the density (and pres¬ 
sure) vanishing at |x| ^ oo rather quickly to guarantee the convergency of all 
integrals involved (whereas the velocity components may even grow). For the 
solution the total mass 


M 



is conserved, the total energy 

= / + dx = £kit) + £p{t) 

R2 
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is nonincreasing function. Let us involve integral functionals 

G{t) = ^J /3|rpdx, 

R2 

Ni{t) = I Xlipdx, i = 1,2, 

]R2 

Ii{t) = I Vipdyi, i = 1,2, 

R2 

Fi{t) = [i = 1,2, 

R2 

where Xi = r = {x,y), X 2 = r_L = ( 7 /, —x). Note that G{t) > 0. 

For the classical (and for the piecewisely smooth) solutions to (El - E3), 


the following relations hold m (see also my- 

G'(t) = Fi(t), (2.1.1) 

Nl{t)=Ii{t),i = 1,2, (2.1.2) 

I[{t) =-Ihit) - phit), (2.1.3) 

l!,it) = lhit)-phit), (2.1.4) 

F[{t) = 2(7 - l)£p{t) + 2£k{t) - lF 2 {t) - pFi{t), (2.1.5) 

F'it) = lFiit)-pF2{t), (2.1.6) 

£'{t) = -2p£k{t). (2.1.7) 


The result can be obtained by means of the Green’ formula. 

Generally speaking, the functions £p{t) and £kit) cannot be expressed 
through G{t), Fi{t), F 2 {t), hit), IVi(t), N 2 {t). But if we choose the 

velocity field with the linear profile V = Ait)r + b(f), we obtain the closed 
system of ODE for finding the coefficients of the matrix A{t) and of the com¬ 
ponents of vector b(t). In general case it is rather complicated and can be 
solved only numerically. 

More simple result can be obtained if 

V = a{t)r -I- /3(f)r_L -|- h{t) = 

= (^ait)I + Pit) 0 ))'■+( 62(0 ) ’ 

with the identity matrix I. It is easy to see that in this case 
hit) = ait)Niit) + Pit)N2it) + biit)M, 
hit) = ait)N2it) - Pit)Niit) + b2it)M, 

Flit) = 2ait)Git) -h biit)Niit) + b 2 it)N 2 it), 

F2it) = 2Pit)Git) + b2it)Niit) - biit)N2it), 

£kit) = ia^it) -h P‘^it))Git) -h ait)biit)Niit) + /3(f)5i(t)iV2(t)+ 
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+a{t)b2{t)N2{t) + l3{t)b2{t)Ni{t) + + bl{t))M, 

£p{t) = KG^-^{t),K = fp( 0 )G^-^( 0 ). 

In that way, all functions involved in the system (2.1.1-2.1.7) are expressed 
through G{t), Ni{t), N 2 {t), a{t), (5{t), bi{t), b 2 {t). 

2.1.2. A{t) of special form, h{t) = 0. The case is simplest. Namely, Fi(t), 
F 2 {t), Ek{f), Ep{t) can be expressed through G{t), a{t), (3{t). Therefore, if we 
denote for the convenience Gi{t) = 1/G{t), then we obtain the system 

G[{t) = -2a{t)Gi{t), 

j3'{t) = a{t){l - 2/3(t)) - ^/3(t), (2.1.9) 

aft) = —a^{t) + 0^{t) — ll3{t) — pLa{f) + (7 — l)KG'l{t). 

For 7 > 1, the functions a{f) and j3{f) are bounded, it follows from the 
expression for the total energy. Really, 

£{t) = (a(t )2 + /3(t)2)G(t) +Tp( 0 )G^-i( 0 )^:^ < £ 1 ( 0 ), 

it implies 

+ < £{0)Gi{t) - £p{0)Gl-"'{0)Gf{t) < + 00 , 

|divV| < 00 . Consequently, the density and pressure are bounded for all 
solutions of the class we want to construct. 

If n = 0, the system (2.1.9) can be integrated and we can obtain 

a{Gi) = ±^C2G'1 - GlGl + {£{{)) - lGi)Gi - F/A, 

m = ciG,{t)+1/2, - / _ = t, 

Jgi{o) 2.Gia{Gi) 

where Cl = C 2 = (a2(0)+C2c2(0)-(£:(0)-/Ci)Gi(0)+/V4)/G7( 0 ). 

In the cases /x > 0 and p, = I = 0, there is a unique stable equilibrium in 
the origin. We have the following asymptotics of the solution components as 
t —)■ 00 : 


if ^ > 0, ^ = 0, thena(t) ~ /3{t) ~ Gi (^ t exp{—pt}, 


27 
Gi(t) 


\2K,^J 


II \ 

^ ^ ^-1/7. 


2 iFi 7 


if ;U > 0,/ / 0, thena(t) ~ -^t /3(t) ~ -^J—t Gi(t) ~ 

27 2 /i 7 


\2Ki,,lj 


-1 


Gi{t)^a\t), m = CiGi{t), 
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where Ki = (7 — l)£p{Q)G\ "'^( 0 ). 

Knowing a{t) and in the standard way we can find components of 

density and entropy from the equations (El), (E3), linear with respect to p 
and p, respectively. We obtain that 

t t t 

|r|,())) = exp (-2 J Q;(r)(ir)/ 5 o(k| exp(- y a(T)(ir), 0 + J j3(T)dT), 


|r|, 0 ) = exp (-27 J Q;(r)(ir)po(kl exp(- y a{T)dT),(j) + J l3{T)dT). 

0 00 

From (E2) and (2.1.9) we get that on the classical solution of the system 
(E1-E3) the relation 

Vp = -{^-l)Gl-^{0)Ep{0)Gj{t)pr 

must be satisfied. Hence it follows that the components of the initial data po 
and po must be asymmetrical and compatible, i.e. connected as follows: 


Vpo = -(7 - l)Gi(0)Ep(0)por- 

For example, one can choose 

1 


Po = 


(l + |r|2)^ 

2a 


a = const > 3, 


1 


Po = 


(7-l)Gi(0)Ep(0)(l + |r|2)«+i- 


( 2 . 1 . 10 ) 

( 2 . 1 . 11 ) 

( 2 . 1 . 12 ) 


2.1.3. A{t) of general form, h{t) = 0. To consider the velocity field (2.1) with 
the matrix 

a{t) b(t) 
c{t) d{t) 


Ait) = 


we introduce the functionals 

1 


1 


1 


Gxit) = T, I Px dx, Gy{t) = - py dx, Gxyit) = - pxydx. 


Then involve the auxiliary variables 

Gi{t) = G^{t)A-<^^+^y\t), G2{t) = Gy{f)A-^^+^^/\f), 

Gs{t) = G^yA-^"'+^'>/\t), (2.1.13) 

where A{t) = Gx{t)Gy{t) — G^yfl) is a positive function on solutions to (El- 
E3). Note that the behaviour of A{t) is governed by the equation 

A'{t) = 2{a{t) + d{t))A{t), (2.1.14) 

and the potential energy Ep{t) is connected with A{f) as follows; 

Ep{t) = Ep(0)AG-i)/2(o)A(-T'+i)G(^)_ 
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To find Gi, G 2 , G3, and the elements of the matrix A{t) we get the system 
of equations 

G'lit) = ((1 -'y)a{t) - {1 + -f)d(t))Gi{t) + 2b{t)G3(t), 

G'2{t) = ((1 - l)d{t) - (1 + 7)a(t))G2(t) + 2c(t)G3(t), 

G'3{t) = c{t)Gi{t) + b{t)G2{t) - + d{t))G3{t), 

a' (t) = —a^{t) — b{t)c{t) + lc{t) — ^.a{t) + K2G2{t), ( 2 . 1 . 15 ) 

b'{t) = -b{t){a{t) + d{t)) + ld{t) - nb{t) - K2G3{t), 
c'{t) = -c{t){a{t) + d{t)) - la{t) - - K2G3{t), 

d'{t) = — b{t)c{t) — lb{t) — iJ,d{t) + K2Gi{t), 

with K 2 = ^(.Bp(0)A(T'-i)/2(0)). 

Now we obtain an analytical result on a behaviour of the coefficients of the 
matrix A{t) as t ^ 00 . 

Proposition 2. 1. Suppose p, = I = 0, that is F = 0. If the system (E1-E3) 
has the solution with the linear profile of veloeity V = ^(t)r, then A{t) ~ |/ 
as t ^ 00 , with the identity matrix I and a constant d > ^. 

Proof. Go to the new variables ai{t) = a{t) — d{t), bi{t) = b{t) + c{t), ci(t) = 
bit) - cit), di(t) = ait) + dit), Gfit) = Gfit) + G 2 it), G^it) = Gfit) - G 2 it). 
In variables ai, 61, ci, di, G3, G4, G5 system (2.1.15) has the form: 


a'fit) = -aiit)diit) - K 2 G 5 it), (2.1.16) 

b'fit) = -6i(t)di(t) - 2K2G3it), (2.1.17) 

c'lit) = -ciit)diit), (2.1.18) 

d'lit) = -]^ia\it) + dlit)) - ]^ib\it) - clit)) + K2G4it), (2.1.19) 

G3it) = —'ydiit)G3it) + -&i(t)G4(t) — -ci(t)G5(t), (2.1.20) 

G'fit) = —7di(t)G4(t) + ai(t)G5(t) + 26i(f)G3(t), (2.1.21) 

G'^it) = -7di(t)G5(t) + aiit))G4it) - 2ci(t)G3(t). (2.1.22) 


Let as t ^ 00 the asymptotics of functions involved in the system be the 
following: ai(t) ~ , 61 (t) ~ L 2 t''^,ciit) ~ L 3 t^^,diit) ~ L4t'"^,G3it) ~ 

Gfit) ~ G^it) ~ M 2 tP^, where Li, i = 1,2,3; Mj, j = 1,2; N are 

certain constants not equal to zero. Note that p 2 < pi and in virtue of A > 0 
the estimate q <pi holds. 

From (2.1.18) we get immediately that 

/3L3f'3-l = -L3L4^'3+'^ 

hence I 4 = —1, L4 = —I 3 . 

Taking the fact into account, from (2.1.16) we have 

- 2K2M2tPL 

The following variants are possible: 
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P 2 < h — 1) hence li = —L4, i.e. li = I 3 , 
P 2 = h — I, hence 


, _ ^ K 2 M 2 _ , 

n — - j — — h — 


From (2.1.17) we have analogously 

l2L2t^^-^ = -L2Lit^^-^ - 


K 2 M 2 

Li 


2 K 2 Nt‘^. 


(2.1.23) 


There are the variants: 

q < I2 — I, hence I2 = —L4, i.e. I2 = I3, 
q = I2 — 1, hence 

, ^ 2K2N , 2K2N 

<2 — —L4 ---— <3---. 

^2 ^2 

From (2.1.19) we get 

-L 4 t -2 = -i(L?t 2 h + L^t- 2 ) _ i(L 2 ^ 2 Z 2 _ ^ 2 ^ 2 Z 3 ) + K2MitP\ 
If li < —1, i = 1, 2,3, then at t —)• 00 


(2.1.24) 


(2.1.25) 


a{t) ~ + Lit ^), 6(f) ~ 

c{t) ~ ^iL 2 t^^ - Lst^^), d{t) ~ ^{Lit~^ - Lit^^), 

that is A{t) ~ Li = const > 0. 

In the case Li = —I 3 > 1, therefore 6 > ^. 

We shall show below that others situations do not appear. 

From (2.1.25) one can deduce the impossibility of the situation I 3 = —1, li < 
— 1, 12 < —l,Pi < —2. Really, in the case I 3 = —Li = —1 — y^l + L3 < —1 or 
I 3 = —Li = —1 — y^l + L3 + 2 K 2 M 1 < —1, it contradicts to the assumption. 

Now let I 3 > —1. Then from (2.1.25) h = I 3 or (and) I 2 = h- Suppose, for 
example, that li = I 3 . Then P 2 < h — 1; from (2.1.22) it follows that 

P2M2tP^~^ = --fLiM2tP^~^ + + 2 L2Nt‘i"^^\ (2.1.26) 


and therefore pi + h = q + h-, h ^ h-, i-e. I 2 > —1- Besides, 


LiMi = - 2 L 2 N. (2.1.27) 

From (2.1.20) we have 

qNti-^ = -jLiNt^-^ + ^L2MitP^+L _ (2.1.28) 

As soon as g — 1 < pi — 1, and pi — 1 < pi + ^2, then pi + ^2 = 7*2 + ^3 < Pi + h, 
and therefore I2 < I3, and as h = I3, then h = I2 = h and pi = q. Besides, 


L 2 M\ — L 3 M 2 - (2.1.29) 

Further, from (2.1.21) we obtain 

piMRPi-i ^ -7L4MifPi-i + - 2L3Atf«+'T (2.1.30) 
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As q + l 3 > Pi-1, then p 2 + li= q + I 3 , p 2 = q, 

L 2 M 1 = 2 L 3 N. (2.1.31) 

From (2.1.27) and (2.1.31) we have ^ and from (2.1.29) we obtain 

^ i.e. L 2 = —F3, L2 = L3 = 0 in spite of the assumption. 

Suppose now that I 2 = h- Then q < I 2 — 1, P 2 < h — From (2.1.20) we 
have q — I < h + Pi, since I 2 > — 1, and therefore ^2 + Pi = ^3 + P2, Pi = 
P2, L2M1 = L3M2. 

Therefore from (2.1.26) we have pi + /i = g + ^2 < Pi + h, h < h, L 2 M 1 = 
- 2 L 2 N. 

From (2.1.30) we get p2 + /i = q + h < Pi + h, P2 = q, h = h, L 1 M 2 = 
—2L^N. In that way, we obtain the contradiction analogous to the previous 
one. 

Now let < —1, li > —1 and (or) I2 > —1. Then from (2.1.25) it follows 
that Pi = 2 li and (or) pi = 2/2 (pi > —2). For example, if pi = 2/i, then from 
(2.1.27) we get p2 - 1 < Pi + ^i, Pi + ^i = q + h < Pi+ ^2, ^1 < h- From (2.1.28) 
P2 + ^3 < Pi - 1, ^2 + Pi > ^1 + Pi > Pi - 1, therefore g - 1 = /2 + Pi- But 
(7 — 1 < Pi — 1, therefore I2 < —1, and ^2 < — 1 in spite of the assumption. 

If Pi = 2/2, then from (2.1.28) we get /3 +P2 < Pi — 1, therefore pi +12 >= 
P2 + ^3, 9 — 1 = Pi + h, h < h- But q — 1 < Pi — 1, therefore I 2 < —1 in spite 
of the supposition. 

So, it remains the unique possibility: < —1, li < —1 and (or) I 2 < —1. 

In the case, as follows from (2.1.25) 

h = -l±yjl- {5iLl + 62 LI - 63 LI - 254 K 2 M 1 ), (2.1.32) 

where 5i = 1, if k = —1, and d* = 0 otherwise, i = 1, 2, 3, (I4 = 1, if p 2 = —2 
and (54 = 0 otherwise. Hence 


5 iLl + 62LI - 53LI - 25 iK 2 Mi < 1, 

(2.1.33) 

if inequality (2.1.33) is strict, then < —1. 

We consider this case, that is < —1, li = —1 and (or) I2 = —1. 

If Zi = —1 / /s, thenp2 = /i —1 = —2. If /2 = —1 / h, then q = I2 — I = —2. 
In that way, in any of these cases we have pi = — 2 . 

At li = —1, I2 < —1, h < —1, from (2.1.26) we obtain p2 = —7T4 + , 

Ij\M\ 

M2 ~ +^^3), 

(2.1.34) 

from (2.1.30) we get pi = —7T4 + 


L1M2 

M. “-P + '-C. 

(2.1.35), 

and from (2.1.24) 

K2M2 , , ^ 

r = t 3 + 1- 

P‘1 

(2.1.36) 
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From (2.1.34), (2.1.36) after excluding Li and M 2 we obtain 

K2M1 = —(2 + '-yh){h + !)• 

As K 2 and Mi are positive, and < —1, then 

h > 

7 

At 7 > 2 it goes already to the contradiction. 

Further, from (2.1.35) 
gether with (2.1.37) gives 


(2.1.37) 


(2.1.38) 


Further, from (2.1.35) and (2.1.36) we have this fact to- 


L? = (2 + 7 / 3 )". (2.1.39) 

From (2.1.32), (2.1.37), (2.1.38) and (2.1.39) we obtain the inequality 


— < -1 - v ' 1-(2 + 7 / 3 ) 2 - 2(2 + 7 / 3)(/3 + 1 ), 
7 


(2.1.40) 


which cannot be true at 7 > 1 . 

If li < —1, I 2 = —1, h < —1, then from (2.1.26) we get p 2 = — 7 A 4 + 

2 L 2 N 


M 2 


= -2 - 7/3, 


from (2.1.28) q = - 7 L 4 + 


L 2 M 1 

2 Ai 


= -2 - 7/3, 


(2.1.41) 


(2.1.42) 


from (2.1.30) 


P2 = -7A4 = 7^3- 


In that way, this fact together with (2.1.41) or (2.1.42) contradicts 

to the fact that L 2 , N and Mi are not equal to zero. 

If li < —1, I 2 = —1, h = —1, then from (2.1.26),(2.1.28) and (2.1.30) we 
get correspondingly 

2L2N LiMi 

+ 4^ = - 2 - 7 / 3 , (2.1.43) 


M 2 M 2 

+2 Ml 


2N 

L 1 M 2 


from (2.1.23), (2.1.24) 


Ml 

772 M 2 2A:2A^ 


= -2 - 7/3, 
= -2 - 7/3, 


= ^3 + 1 - 


(2.1.44) 

(2.1.45) 

(2.1.46) 


Li L 2 

In that way, multiplying (2.1.43) by (2.1.45), taking into account (2.1.46) we 
get 

M 2 LI 


(2 + 7/3) — Li + 


Ml 
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from (2.1.44), (2.1.45), (2.1.46) 

Ml = Ml, (2.1.45) 

that is 

{2 + ^hf = Ll±Ll (2.1.48) 

As above, from (2.1.34), (2.1.37), (2.1.38) and (2.1.48) we get the inequality 
(2.1.40) which cannot hold at 7 > 1. 

It remains the unique possibility = ^2 = ^3 = — 1- In the case p 2 < 
—2, q < —2 and, as follows from (2.1.26), pi = p 2 or pi = q, that is pi < —2. 
Therefore, as follows from (2.1.32) 

lI + lI-LI = 1. (2.1.49) 

If Pi = p 2 , then IM 2 I < Ml. If pi = P 2 , q < PI 7 then from (2.1.26), (2.1.30) 
we obtain that P 2 + 7 = Pi + 7 = hence Ml = M|, it goes to the 

contradiction. 

If Pi = q, P 2 < Pi, then from (2.1.26), (2.1.28), (2.1.30) we get 

LiMi = - 2 L 2 N, A = (2.1.50) 

where A = g + 7 = pi + 7 . It follows, in particular, that = — ALi, L 1 L 3 = 
— AL 2 , A^ = —L 2 T 3 . Besides, if Gi{t) ~ Nif^^, G 2 {t) ~ N 2 t‘^^, t 00 , where 
Ni,N 2 are some positive constants, qi = q 2 = Pi = Q, A^i = ^ 2 , Mi = 2A'i, 
and Nl > N'^ in virtue of Gi{t)G 2 {t) — G^t) > 0. In such way, from (2.1.50) 
we get that > A^, A^ > and taking into account (2.1.49), 

< 1, < 1, A^ < 1. That is if A > 1 (7 > 3), the conditions mentioned in 

the paragraph cannot hold together. 

At last, consider the casepi = P 2 = Q- Then from (2.1.26), (2.1.28), (2.1.30) 
we have A = that is the system 

of linear homogeneous equations with respect to the variables Mi, M 2 , N 

LiMi — XM 2 + 2 L 2 Xf ~ 0 ) 

£ 12^1 — L 3 M 2 — 2XN = 0, 

XMi — L 1 M 2 2 L 3 N = 0. 

For the existence of its nontrivial solution the determinant of the system must 
be equal to zero, i.e. 

ATi — 2AT2A3 + L1L2 + T1T3 — A^ = 0. 

Taking into account (2.1.49), involve the function with respect to the variables 
L 2 and L 3 , where A plays the role of parameter, namely 

^\{L 2 , A3) = A (1 — lI + A3) — 2AA2A3 + (L2 + A|) Y^l — L2 + Ag — A^. 

By the standard methods one can show that the function is not equal to zero 
at A > 1, i.e. at 7 > 3 (remember, that pi < — 2 ). 

So, it remains to investigate the case 1 < 7 < 3. 
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Consider equation (2.1.14), which can be written as 

A'{t) = 2di{t)A{t). 

If we suppose that A{t) ~ const ■t'^, t —)• oo, m is a constant, then m = 2 L 4 = 

2. From (2.1.14) we can get Ep{t) ~ const ■ 0), t —)• 00 . Therefore, if 

we denote E the quantity of the total energy of the system, then Ek{t) = E — 
Ep(t) E{1 — Cit^ '•'), here and further Ci are some positive constants. But 

Ek{t) ~ C 2 Gk{t)t~‘^A'^ ~ where G^ is at least one of functions 

Gi, G 2 or G 3 . That is G^ ~ G 4 (l — Git^~"<)t^~'^ ~ G^G~'^. But then 1 — 7 < 

Pi < —2, 7 > 3, and 7 does not belong to the interval under consideration. 

So the proof of the proposition is over. 

Remark 2.1.1 One can show more shortly, that in the physical case 1 < 

7 < 2 the situation Zi = ^2 = ^3 = if impossible. Taking into account 
(2.1.13) we have G 1 G 2 — G| = A~'^ ~ const ■ t~‘^'^, t —)■ 00 . But the degree of 
the leading term of the expression G 1 G 2 — G 3 is not greater then 2pi, therefore 
Pi > — 7 , as Pi < —2, then 7 > 2. 

Remark 2.1.2 Actually I3 = — L 4 = —2, G 1 G 2 —G 3 ~ const ■ t~‘^'^ , t 00 . 

2.1.4. A{t) of special form, b(t) / 0. Suppose that the velocity field is of form 
( 2 . 1 . 8 ). 

In the case we need to analyze rather complicated system of 7 equations: 
G'{t) = 2a{t)G{t) + bi{t)Ni{t) + b 2 {t)N 2 it), 

N[{t) = a{t)Ni{t) + l 3 {t)N 2 {t) + bi{t)M, 

N2{t) = a{t)N2{t) - l3{t)Ni{t) + b2{t)M, 

{a{t)Ni{t) + l 3 {t)N 2 {t) + bi{t)My = 

-p{a{t)Ni{t) + l 3 {t)N 2 {t) + bi{t)M) - l{a{t)N 2 {t) - (3{t)Ni{t) + b 2 {t)M), 
{a{t)N 2 {t) - mNi{t) + b 2 {t)M)' = 

-p{a{f)N2{t) - f3{t)Ni{t) + b2{t)M) + l{a{t)Ni{t) + (3{t)N2{t) + bi(t)M), 
(bi(t)Ni(t) + b 2 (t)N 2 (t) + 2a(t)G(t)y = 

-fj,( 2 a(t)G(t)+bi(t)Ni(t)+b 2 (t)N 2 (t))+l( 2 / 3 (t)G(t)+b 2 (t)Ni(t)-bi(t)N 2 (t))+ 
(blit) + bl{t))M + 2(7 - l)KG^-^{t) + 2 G{t){a\t) + ^\t))+ 
2a{t)bi{t)Ni{t) + 2 f 3 {t)bi{t)N 2 {t) + 2 a{f)b 2 {t)N 2 {t) - 2 j 5 {t)b 2 {t)N 2 {t), 

{b 2 {t)Ni{t) - bi{t)N 2 {t) + 2/3(t)G(t))' = 
-l{ 2 a{t)G{t)+bi{t)Ni{t)+b 2 {t)N 2 {t))-p{ 2 l 3 {f)G{f)+b 2 {t)Ni{t)-bi{f)N 2 {t)). 



SMOOTH SOLUTIONS TO THE MULTIDIMENSIONAL BALANCE LAWS 13 

We could write the normal form of the system, but it is very long. Mention 
only that we can always express the derivatives of a, (3, bi, 62 explicitly, because 
the determinant of the corresponding algebraic system is equal to 

As the Holder inequality shows, it is positive for the solutions we consider. 
The compatibility condition in the case has the form 

Vpo = Po(cir + C2r_L + Cq) 
with scalar constants ci, C 2 and constant vector cq. 


2.1.5. A{t) of general form, b(t) 7 ^ 0. The corresponding system consists of 
11 equations for unknown functions a{t), b{t), c{t), d{t), bi{t), b 2 {t), Gx{t), 
Gy{t), Gxy{t), Ni{t), N2{t). It seems that it can be analyzed only numerically. 

As soon as the components of velocity field are found, we can always find 
the density and the pressure. Namely, 


r r 

p{t,x) = exp{— J ii A{t) dr) pq{M~^ {t)x — j M“^(r)b(T)dT), (2.1.51) 

0 0 

t t 

p(t, x) = exp (—7 J trA{T)dT)po{M~^{t)x — J M“^(r)b(r)dr), (2.1.52) 

0 0 

where M{t) is a solution to the matrix equation = A{t)X{t) |18j . 

The compatibility condition in the case has the form 


Vpo(x) = po{x){Cx + Co) 


with a constant matrix C and a constant vector cq. 


2.2. The system in the physical space (n = 3). Now we consider the 
right-hand side of (El), corresponding to the Coriolis force and the Rayleigh 
friction in the real physical space. Namely, 

/ 0 Ws —U2 \ 

f =-^V + (5[V X w] = (-/rl + 5 -ws 0 wi )V, (2.2.1) 

y UJ 2 -oji 0 j 

where d = 0 or 1 , a; is a constant vector (wi,ti; 2 ,^ 3 ), /r > 0 is a constant, 
[a X b] is the cross product of the vectors a and b. 

Besides of G{t) = 5 / p|rp dx and Fi{t) = f (V,r)pdx, the generalization 
K3 R3 

of the functionals, involved in 2.1.1, to the 3D case, we consider 

m) = [ 

R3 


(V, [uj X r])pdx, 
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F'sit) = /([V X cj], [oj X r])pdyL, 

IR3 

H{t) = ^ J ([a; X r])Vdx. 

R3 

For the smooth solutions to the system (E1-E3) with the right-hand side 
(2.2.1) we have 

G'(t) = Fi(t), 

F[{t) = 2£k{t) + 3(7 - l)£p{t) + 5F2{t) - /uFi(t), 

F’{t) = SFsit) - pF2{t), 

H'{t) = 0. 

We seek the solution with a special linear velocity profile 

/ 0 Ws -012 \ 

V = Q;(t)r -|- /3{t)[r x w] = {a{t)E -|- j3{t) j —ws 0 oji 1 )r. 

\ UJ2 —UJl 0 J 

In this case 

Fi{t) = 2a{t)G{t), 

F2it) = -2f3{t)H{0), 

Fsit) = —2aF[{0), 

£k{t) = a^it)G{t) + p\t)H{0), 

3(7-1) 

£p{t) = const ■ G 2 (t). 

Changing G~^{t) = Gi{t), we obtain the closed system of ODE 
a'{t) = -a‘^{t)-i2a{t)+/3^{t)Gi{t)H{0)-S/3{t)Gi{t)H{0)+ ^^^ ~ KG^^{t), 
/3'{t) = 5a{t) — pPit), G'i{t) = —2a{t)Gi{t), 

^ , s ^.3(7-1) 

with K = £p{0)G 2 (0). Since the total energy is nonincreasing, we can 

draw similarly to the 2D case that |a(t)| < 00, \f3{t)\ < 00. 

In the situation with p = 5 = 0, with p > 0, 5 = 0 and with p > 0, 6 = 1, 
there is a stable equilibrium in the origin a{t) = I3{t) = Gi{t) = 0. 

If = 0, 5 = 0, we have 


a(t) ~ t , P{t) = P{0) = const, Gi{t) ^ const ■ t ^,t^oo. 

If ^ > 0, (5 = 0, then 

1 1 2 

~ -7^ , P{t) =/3{0) ex.p{—pt),Gi{t) ^ const ■ t 37 - 17 ^ 00 . 

37 — 1 

In the case ^ > 0, <5 = 1, 

1 


a{t) 


37-1 


t ,/3{t) 


1 


^(37 - 1 ) 


t ,Gi{t) ^ const ■ t 37-i,t^oo. 
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What about other component of solution, they can be found by the same 
formulas, as in the 2D case ((2.1.50 - 2.1.51)), and must satisfy the compati¬ 
bility condition 

Vpo = -(7 - l)G'i(0)£;p(0)por, 

the same as (2.1.10). 

Remark 2.2.1 We can also consider the velocity with linear profile of 
general form, but the system of integral functionals is complicated and it is 
difficult to analyze it. 

2.3. Theorem on the interior solutions. Further we need to obtain the 
symmetric form of system (E1-E3). Eor this purpose we consider the entropy 
S, connected with the components of solution by the state equation p = . 

Thus, instead of (E1-E3) we obtain other system 

-I- div (pV) = 0, (2.3.1) 

pdtV + {pY, V) V + Vp = pf (x, t, V, p, S), (2.3.2) 

dtS+{Y,YS) = 0. (2.3.3) 

Here /?(t,x), V(t,x), S(t,x) are components of the solution, given in M+ x 
M"', n > 1 (density, velocity and entropy, respectively). 

The systems (E1-E3) and (2.3.1-2.3.3) are equivalent for p > 0. 

Put the Cauchy problem for (2.3.1-2.3.3): 

p(0,x) = po(x) > 0, V(0,x) = Vo(x), 5(0,x) = 5'o(x). (2.3.4) 

Denote po(x) = e'^°^^^/?Q(x), and remark that the first of the Cauchy con¬ 
ditions can be replaced to p(0,x) = po{Y) > 0. 


Definition. We shall call the global-in-time classical solution (p(t,x), 
V(t,x), S{t,x.)) to the system (2.3.1-2.3.3) the interior solution, if 

(p^, V, S) E C^(M+ X M’^), 

where p = e^p^, and any solution (/9(t,x), V(t,x), S'(t,x)) to the Cauchy 
problem (2.3.1-2.3.4) with the sufficiently small norm 

-p^'^"^^/^^(0,x), Vo(x) - V(0,x), 5'o(x) - 5(0,x)||j^m(Kn), 

7— 1 

m > n/2 -|- 1, is also smooth, such that {p 27 , V, S') E C'^(M+ x R"’). 


Note that the trivial solution is not interior at least for F = f = 0. But 
a set of the interior solutions is not empty. In |15j (a generalization of m) 
for f = 0 it was shown that the solution (0, V(f,x), const) to (2.3.1-2.3.4) 
is interior, provided V(t,x) is a globally smooth solution to the equation 
dtY + (V,V)V = 0 such that the spectrum of its Jacobian is separated 
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initially from the real negative semi-axis, -DV(0,x) G Ii)^V(0, x) G 

denote by the vector of all spatial derivatives of the order 

k). 

The result is clear from a physical point of view: the velocity field with a 
positive divergency spreads the initially concentrated small mass, that prevents 
the singularity formation. 

As we show below, some of solutions with linear profile of velocity are 
interior with the density not close to zero. 


For the solutions to (2.3.1--2.3.3) with a finite moment we have inf p = 
infp = 0, thus we need to use the symmetrization proposed in [T]. After 

7-1 

involving the new variable H = up , k = —-—, we obtain the symmetric 

7-1 

form of the system (2.3.1-2.3.3) with the solution U = (II, V, S): 

S 1 s 

+ (V, V))n + l^exp(-)ndiu V = 0, (2.3.5) 

7 27 


{dt + (V, v))v + ^ exp(-)nvn = fi(t, x, n, v, 5 ), (2.3.6) 

2 7 

{dt + {\,V))S = Q, (2.3.7) 

where fi = f(t,x,e-f (3)7^ ,V, 5). Let Hq := np}^, Uq := (Hq, Vo,5o). 
In this way, the Cauchy problem (2.3.1-2.3.4) is transformed to the problem 
(2.3.5-2.3.7), with the initial data U(0,x) = Uo(x). Further, suppose that 
the system (2.3.5-2.3. 7) has a classical solution U := (n, V, S), where V := 

- 7—1 - 

A{t)r -I- b(t), n := Kp 27 . Denote also U — U := u := (tt, v, s). 

Before formulating Theorem 2.1 we perform certain transformations. Firstly 
we go on to the system with the solution u. Then, following m, we carry out 
the nondegenerate change of variables such that the infinite semi-axis of time 
turns to semi-interval. 

Assume that there exists a nondegenerate square (n x n) matrix Ai{t) 
such that Ai(t)(A)"^(t))' = A{t), moreover, Ai{t)A{t) = A{t)Ai{t). Choose a 
positive decreasing function A(t) such that the integral \{T)dT converges 
to the finite value cToo, and set a{t) = Jq A(r)cir. Let (cr, y) := {a{t), Ai(t)x.) 
be new variables. Then Vx = ^fVy, diVxV = divyAiV, = X~^(t){dt + 
^rVx). In that way, the semi-infinite axis of time goes to the semi-interval 
[0,(Too)- Note that there exists the inverse function t = t{a). Further, involve 
the variables W = X~^{t)Ai{t)v, P = A“'?(t)7r, the constant q will be defined 
below. Let U{a, y) := (P, W, s), P := A-^fi, W := A-^AiV, U := (P, W, S). 
So we get a system 


(a. + (W,Vy))P + 


7-1 


(P + P)divy W = 


2 
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-(W, Vy)P - PQi - X-^A^h, VyP), (2.3.8) 

(d^ + (W, Vy)) W + a){P + P)VyP = 

a){P + P)VyP + ^^( 5 , a)PVyP- 

-QaW + G- X-\Aih, Vy)W, (2.3.9) 

{d^ + (W, Vy))s = -(W, Vy)5 - A-^^ib, Vys). (2.3.10) 

Here we take into account that t = t{a), x = A^^{t{a))y, A = X{t{a)),A = 
A{t{a)),Ai = Hi(t((T)),b = b(t(cr)) and denote 

Qi = Qi{t[a)) := A(t)"^(^^^y^trH(t) + g(ln A(t))'), 

Q2 = Q2{t{a)) := X{t)-\{\nX{t))'E + A{t) + Ai{t)A{t)A-{\t)), 

'k(5, cj) := exp( —)i?H 2 , 

7 

= AMa)) := Hi(t)Ht(t)(detHi(t))-2/- 

R = R{t{a)) := X^^-'^{t){dei Ai{t)fl^, 

G = G{a, y,U,U)-=X-^ {t)Ai {t) (fi {t, x, A"? (t) {P + P), 
X{t)A^\t){W + W),{S + s)) - fi(t,x, A^(t)P, X{t)A^\t)W, S)). 

Note that A2 is a family of invertible matrices. 

The initial data Uq and the vector-function U induce the Cauchy data 
V(o{y) for the system (2.3.8-2.3.10) as follows: Uo{y) = A“'^(no(A)"^y) — 

n(0,H-iy), A-iHi(Vo(H-V)- V(0,H-V)), SoiA'^y) - S{0, A-^y))\^=o- 

After multiplying (2.3.9) by 'I'“^(S', fi) the system becomes symmetric hy¬ 
perbolic. We can apply the theorem on a local existence of the unique so¬ 
lution to the Cauchy problem for symmetric hyperbolic systems [3] to the 
problem (2.3.8-2.3.10), Z^(0,y) = Uo{y), provided 'I'“^(S',(t) is uniformly 
positive. Sometimes T“^(S',a) is really uniformly positive (see [M]), and 
in that cases we can conclude that if Uq € H^(MP),m > 1 -|- n/2, then 
U G nf^QG^{[0,a^y,H'^-^{R^)),a^ > 0. Since a* ^ 00 , as \\Uo\\h^ ^ 0, 
then choosing the initial data Uq small in the Sobolev H^— norm, we can 
extend the time of existence of the smooth solution U to the Cauchy problem 
up to fj* > CToo- In this case we conclude that U G nJhQC''^([0, (Too); H'^~^{R^)) 
and u G n^oCT([0,oo);F”*-^'(M’")). 

However, generally speaking, T“^(S', cr) —0 as t — 00 , therefore we can 
extend the time of existence of the smooth solution U to the Cauchy problem 
up to [0,(T*), for any cr* < (Too, choosing the initial data small in the Sobolev 
H^— norm, but we cannot guarantee the existence of the solution for [0, (Too)- 
Thus, we need to act differently. 

If t* is the supremum of the time of existence of the smooth solution U for 
the Cauchy problem (2.3.5-2.3.7), U(0,x) = Uo(x) (the solution u preserves 
the smoothness during the same time), then If < (Too{t^ < 






18 


OLGA S. ROZANOVA 


oo), then lim sup (||Z^||l^(cj) + \\^yU\\L^{a)) = +oo (lim sup(|| u||l^( t) + 

(T — t — 

llVxulli,^ (t)) = oo). Therefore, if we prove that ||Z^||/fm(cj) < +oo,(T G 
[0,(Too), then we establish that ct* = aoo{t* = oo). 

Let 



Further, fiv stands for the matrix ||^^||, Vnfi and V^fi for the derivatives 
of fi with respect to 11 and S, correspondingly. 


Theorem 2. 1. Let the function fi(t, x, II, V, 5) have the derivatives with 
respect to all arguments up to the order m + l,m > n/2 + 1. Suppose that 
the system (2.3.1-2.3.3) has a global in time classical solution (p, V,,?) with 
linear profile of velocity V = A(f)r + b(t) such that 

a) p^{t,iL) G nf+o^C'^([0,oo);F"^-^'(M")); 

DS{L^) G n”LoCJ'([0,oo);iL™-^'(M’^)); 

b) there exists a matrix Ai{t) such that A{t) = Ai{t){Af^(t))', A{t)Ai{t) = 
Ai{t)A{t), detail(t) > 0, fort >tQ>0. 

Let us assume that there exist a smooth real-valued decreasing nonnegative 
function X{t), a constant q and a matrix U(j){t) with real coefficients, having 
the following properties: 



\{T)dT < oo. 


(2.3.11) 



A'^(r)(det yli(r))^'^'^ dr < 


oo. 


(2.3.12) 


R{t)' > 0 and Qi{t) > 0 for t >to >0, (2.3.13) 


A 2 ^{t)Uff,{t) is a skew-symmetric matrix, 

the following functions, vector-functions and matrices are bounded in [to^oo) x 
(under bounded (11, V, 5)); 

A“'^(t) exp(—-tr^(t)), exp(— f A{T)dT), (2.3.14) 

2 Jto 


A ^{t)h{t), A ^{t)A2^{t)A2t) (2.3.15) 

A-^Zl-^(Vnfi), A-(''+i)zi-^(Vsfi), (2.3.16) 

Qs ■■= {Q2{t) - X-\t)D'^{Ai{t)Vvii)Af\t) - U^{t)), \ J\ = 0,1,... ,m. 
Then the solution {p, V, S) is interior. Moreover, 


(p(7-i)/27 _^(7 -i)/27^ V - V, 5 - 5) G n™ o<^^'([0,oo);/7™-^’(M^)). 



SMOOTH SOLUTIONS TO THE MULTIDIMENSIONAL BALANCE LAWS 


19 


Remark 2.3.1. If Qs depends only on time, then we can require only the 
nonnegativity of this expression for t > to. 

Proof. Let a G [0,cr*). Suppose for the sake of simplicity that to = 0, 
otherwise we can consider the Cauchy problem at ao = cr(to) < cr* with the 
initial data induced by the Cauchy data Uq. Note that choosing Uo sufficiently 
small in the - norm we can prolong the smooth solution U up to do. Denote 

{X, Y, Z) : — )• R X R" x R a vector-function from Introduce a 

norm 

[X, Y, Zf{a) := J{{X^ + Z^) + a)Y)dy. 

R" 

It is equivalent to the usual - norm, but not uniformly in time. 

For any p G Y we define 

- m 

W = 2 E D^s]Ha), FM = 

\J\=p p=o 

Let us compute 

^ = X] [/ {D^PD-^d^Pdy + j D-^sD-^d^sdy+ 

DI=pr" R" 

+ j D'^W'if-\S,a)D-^daWdy+ 

R" 

+ j- i^-i(5,d)5,5)D^Wdy] = 

R" 

= II -I- /2 + Is + h, 

the integrals 1^, k = 1,2,3,4, are numbered in the consecutive order. 

We estimate every of the integrals in the standard way using the Holder and 
the Galiardo-Nirinberg inequalities (see [IB],[TS]) under assumption (2.3.11 - 
2.3.16) of Theorem 1. We denote by c* certain positive constants, which do 
not depend on a. 

Let us begin from R. 

R = - y D“RD“((W,Vy)R+^^(R+R)divyW+(W,Vy)R+PQi(Rd)))dy. 

Rn 

The integral from the terms of higher order D°‘PD°‘{{'W, Xy)P+I—^ —RdivyW) 
can be reduced by integration by parts to 
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7-1 


J D"P{(VyP,Il“W) + (PdivyZl"W - D"(PdivyW))}dy := In. 


Rn 


One can show applying the Holder and the Galiardo-Nirinberg inequalities 
that 

Ill < CillW||2(||W||2||VyW||oo + ||P^W||2||VyP||oo) < 

(the details in [l6],[T5].) Further, according to (2.3.5), (2.3.13) and (2.3.14) 


Zl"PP)"(WVyP)dy 


\Rn 


< ‘'exp(-^^^y^trH)Fm < C2P^/^Pm, 


P)“PZl“(P)divyW)(iy 


Rn 




D^PD^{PQi)dy > 0. 


Rn 


The integral R can be estimated analogously: 

h < c^R^/'^FU^Ep + c^R^/^F^. 

Further, 

Is = y (P“W)*T-^(P, a) {(W, Vy)P“W - P“((W, VyW)} dy+ 

Rn 


+ 


i j |(P“W)*T-i(5,cj)P“W)divyW- ■^^^^y^(P"W)*T-^(5, (t)P"W j dyT 


4 


+ 


Rn 
7-1 


{ {D^W)*^-\S, a)D^ (PT(S, a)VyP + T(5, a)Vy(PP)) } dy+ 


Rn 


+ ^ J {(P“W)*T-i(5,a)P“((T(5,a)-T(S,a))Vy(p2))}ly+ 

Rn 

- J (P“W)*T-i(5, a)P“ (Q 2 W - G) ly. 

Rn 

First two integrals can be estimated from above by the values 

C6P-'||PPW||i||VW||oo + 

C7P-'||PPW||i||W||2 < c^R^/^FU^Ep, 
the third, containing the terms of the form 

k k 

DPW{YI{D^Sfi)D^-’^PDP+^-^P, DPW{YI{D^Sfi)D^-'^PDP+^-^P, 
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or 


k 

DPW{Yl{D^Sf^)D^-^PDP+^-^P, l<l<p,l<k<l,'^jPj =k, 
j=i j 

as follows from the Galiardo-Nirinberg inequality and (2.3.14), can be esti¬ 
mated by the sum of form 


>0. 

i=9 

The fourth integral contains the term of form 
k 

exp(0s )DVW([\{D^SfpD^-^PDP+^-^P, l<l<p,l<k<l, 
i=i 


= fc,|0|<l, 

j 

and can be estimated in a similar way, taking into account that the value of 
||s||oo is bounded. 

The last integral taking into account condition (2.3.16) can be estimated 
by the quantity cgFm- Remark here that 

G = X-\t)Ai (t)Vvfi(t, X, n + 01TT, V + 02V, S + 03 s)A^^ (t)W+ 
A«-2(t)Vnfi (t, X, n + 0i7r, V + 02V, S + 03s)R+ 

A-2(t)Vnfi (t, X, n + 0i7r, V + 02V, S + 03s)s, 
where 0* G (0,1), i = 1, 2,3, and 

J {D^Wy^-\S,a)D^{U4t)W)dy = 0 , 


since 

(w,y42^(t)t/<j(,(t)w) = 0 
for any real-valued vector w. 

Further we estimate I 4 . Remark first of all, that according to (2.3.13) 
R'{a)/R{a) > 0. Taking into account (2.3.15) we obtain 


9T-i(5,cj) 

Thus, 


T-2(5, a)T(5, a)R' 
R 


+ 0(T-n5,a)) <ci2T-n5,a). 


h < C13(1 + llWVySiloo + ||WVy5||oo)^p < Ci3(l + R^^^Fm + R^/^F^^)Ep. 


So we get 

h < ci{Ep + RP^EP^Ep + RP^Fm), 
h < C2{R^/^FU^Ep + RP^F^), 
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4 


h < c^{Fm + + R^/^Fm + R^/^eF^ ^ 7 * > 0 , 


2=1 


h < C5(l + RF^Em + RF^EF^)Ep. 



Fm 7 CQ{Fm + R^F(^P^ -|- E Fm ' ) ) 7 * — 


j=l 


Set Km{cr) = e Em- Then 


6 


A'™ < C7(A^ + ^(A™)3/2+7)iji/2^^^ > 0, (2.3.17) 


2 = 1 

where the constant cy depends on cg, cToo- Let 

9 



5 > 0. The integral diverges in the zero, such that 0(0) = —oo. Integrating 
inequality (2.3.17) over a we obtain 


cr 



0 


moreover, as follows from (2.3.12), the integral in the right-hand side of the 
last inequality converges as u —>■ cToo to the constant C, depending only on 
the initial data. Choosing Am(0) (that is the Hm— norm of the initial data) 
sufficiently small, one can make the value 0(Am(O)) -|- C later then 0(-|-oo), 
it signifies that Am(o') and Fm{o') are bounded from above for all a G [0, fioo) 
and (T* = (Too- So, u G nJLQC-^([0, oo); (M"’)) and Theorem 2.1 is proved. 

Remark 2.3.2 The solutions with linear profile of velocity, satisfying the 
conditions of Theorem 2.1, exist. In the case investigated in [E], f = 0, 7 < 


1 + ^, A{t) = {E + tA{0))-^A{0), SpA(0)nM_ = 0,b(t) =0, Ai{t) = A{t), 
p = 0, S = const. Here A = (1 -|- t)“^, q = , 


2.3.1. Corollary and examples. The following Corollary helps us to prove that 
among the solutions to the Euler system constructed in section 2.1 and 2.2 
there are interior ones. 

Corollary 2. 1. Let? = L{t)\', where L{t) = —pi-\-Ui{t), p is a nonnegative 
constant, Ui (t) is a one-parameter family of skew-symmetric matrices with the 
smooth coefficients. Assume that the system (2.3.1-2.3.3) has a global in time 
classical solution {p, A{t)r, S), satisfying the condition a) of Theorem 1 and 
A{t) ~ a{t)I -\- (d{t)U 2 , as t ^ oo, where a.{t) and I3{t) are the real-valued 
functions, U 2 is a skew-symmetric matrix with constant coefficients. Moreover, 
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suppose that a{t) = with d = const, such that <5 > 0 for /U > 0, and S > ^ 
for n = 0. Then this solution is interior. 

Proof. Here h{t) = 0. The matrix Ai{f) = exp(—yl(r)(ir) satisfying 

the condition b) of Theorem 2.1 exists, detHi(t) = exp(—n a{T)dT), as the 
eigenvalues of a skew-symmetric real-valued matrix have the real part equal to 
zero. The first and second conditions in (2.3.16) hold, since f is independent 
on p and S, Q 3 depends only on t{see Remark 2.3.1). We choose as U^{t) 
the matrix A 2 {t)Ai{t){j 3 {t)U 2 — Ui{t))Af^{t) (such that Af^{t)U^{t) is skew- 
symmetric) to obtain Qs = A“^(t)((ln A(t))' -|- 2a{t) + pA 2 {t))I. If \{t) = 
exp(—e J^^a{T)dT),e = const > 0, £<5 > 1, then both integrals in (2.3.11), 
(2.3.12) converge. To guarantee the nonnegativity of R'{t),Qi{t),Q 3 {t) we 
must satisfy the following inequalities for e and q: q < 1, e > < 

2 q^'^ ’ £ < 2 (the last inequality arises only in the case p = 0), the conditions 
(2.3.12) hold for e < /U > 0, it implies , y < 

a < q < q := min{^^lif^, 2 +i]~li) }, and we choose q G {0,q). If 

p = 0, then additionally we have inequalities ^ < £ < 2, < £ < 2. Since 

<5 < i, then we can choose £ to satisfy the first inequality, for the second we 
have g < i. In this case we choose q G (0, min{g, i}). Thus, all conditions of 
Theorem 2.1 are satisfied and the Corollary is proved. 

Examples. Thus, as follows from Corollary 2.1, in the cases p = I = 0 
and /i > 0 in the Sections 2.1.1 (n = 2) and 2.2 (n = 3) we have constructed 
the velocity field for the interior solution. Really, the right-hand side has 
the form indicated in the Corollary statement, the velocity has linear profile 
V = A{f)r, with A{f) of special form, that is A{f) = a{f)I + /3{t)U2 with the 
skew-symmetric U 2 . For p = I = 0, we have a{t) ~ t~^,t 00, 5 = 1 > ^. 

For p > 0, we get a{t) ~ ^ oo,d = ^ > 0 {n = 2), and S = 

3 ^^ > 0 (n = 3) (remark that the denotation 5 has here other sense than in 
the subsection 2 . 2 .) 

The solution constructed in the Section 2.1.2 is also interior for p = I = 0. 
As follows from Proposition 2 . 1 , in the case A{t) ~ a{t)I, where a{t) = 
5t~^,t —>• 00 , (5 > 5 . 

Moreover, as the numerical analysis suggests, all solutions with linear profile 
of velocity (of general form) are interior if /x > 0 , Z = 0 , either for b(t) = 0 or 
b(t) 7 ^ 0 (it seems that in the last case h{t) ~ with the constant vector 

B). However, the last assertion can be considered only as a hypothesis. The 
case /i > 0 , Z / 0 is more complicated, because it seems that it accepts stable 
equilibriums except of the origin. 

Remark 2.3.3 We can find the entropy function for the solution con¬ 
structed in Sections 2.1 and 2.2. For example, in the simplest case (Section 
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2.2.1) 

t t 

5'(t, |rl,(/>) =S'o(|r|exp(-y a(r)(ir), (/> + J l3{T)dT). 

0 0 

For example, for the initial data (2.1.12) Sq = const + (a (7 —1)+7) ln(l + |rp). 

In spite of the components of the pressure and the density must vanish as 
|x| —)• oo, the entropy may even increase (remember that the conditions to 
Theorem 2.1 require only the boundedness of the entropy gradient). 

Remark 2.3.4 In the case G'(O) ^ 0 the density cannot be compactly 

27 2 _ S 

supported (in contrast with G(0) = 0). Really, since p = , p = Tn'-i e r , 

_ s 

then due to the compatibility conditions ttVtt ~ const ■ e t" , |x| —)■ c — 0, 
where c is a point of the support of tt. Therefore, for C^— smooth tt it occurs 
that S —)• + 00 , |x| c — 0, and we cannot choose any smooth initial data. 
It is interesting that if one requires only C^— smoothness of tt (vr ~ const ■ 
(c — |x l)‘'^ |x| —)■ c — 0, TT = 0, |x| > c > —0), the condition may be fulfilled. 
Moreover, p and p will be of the C^- class of smoothness, however, neither the 
theorem on the local in time existence of the smooth solution, no Theorem 2.1 
can be applied. 

3. Acceptable velocities: generalization of velocity with linear 

PROFILE 

3.1. The auxiliary system of transport equations. Consider the follow¬ 
ing system of equations associated with (E1-E3): 

atV + (V,V)V = 0. (3.1.1) 

Here V is a vector field from the tangent bundle of S. 

We look for the smooth solution to (3.1.1) of the separated form 

V = A(t)A(x) (3.1.2) 

with a square (n x n) matrix A{t). 

The following possibilities can be realized: 

I. 

VjA-^ = const ■ Sj, (3.1.3)(Al) 

with the Kronekker symbol 6j, A'it) = — A^(t). 

If Sp{A{0)) n M_ = 0, then A{t) is bounded all over the time t>oo. 

It is easy to see that for the euclidean metrics there are the unique possi¬ 
bility A = or -|- 6, where r is the radius-vector of point, a, b are constants, so 
we come back to the velocity with the linear profile. 

II. 

A* = A^VjA\ i = l,...,n, (3.1.4)(A2) 

and A{t) = a{t)I, with the identity matrix I, a{t) is a function such that 
a'{t) = —a?{t). 
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In the second case we essentially restrict the set of matrices, but the set of 
corresponding vector fields A is more rich then in the case (Al). 

Remark 3.1.1. The equations (Al) may be incompatible (f.e. for the 2D 
sphere with natural coordinates). 

Remark 3.1.2. If A satisfies (Al), then it satisfies (A2). 

3.1.1. The vector field satisfying condition (A2) with a constant divergency. In 
our further considerations the vector fields A satisfying (A2) having a constant 
divergency, play an important role, therefore we study them in detail. 

Theorem 3. 1. If the nontrivial potential vector field satisfying (A2) has a 
eonstant divergeney D = Dq, then Dq > 0. 

Proof. Really, 

(VA) = Vi(A^ VjA*) = 2 ^ ViA^ViA^' + ^(ViA*)^ + ^ A^Vj(VA). 
id i¥=j i j 

(3.1.5) 

(in the formula we do not sum over indices i in (VjA*)). 

If the field A is potential (that is there exist a function <I> such that A = 
V^h), then VjA-^ = VjA*. If we suppose that the divergency is constant, then 

Do = 2^(ViA^')2 + ^(V,A*)2. 

i¥=j i 

Thus, Dq is a sum of squares and it is therefore nonnegative. However, if 
D = VA = 0, then VjA-^ = 0 under all combinations of indices, therefore A 
is constant and, as follows from (A2), equal to zero, that is trivial in spite of 
the assumption. 

Theorem 3.1 has the evident corollary. 

Corollary 3. 1. There is no nontrivial potential divergency free veetor field 
satisfying (A2). 

Remark 3.1.3. As we shall see below, the possibility to choose the diver¬ 
gency free vector field should involve the possibility to construct the solutions 
to the Euler equations with the properties of localization of mass in a point. 
But the situation is restricted for the systems, containing pressure. 

The case of dimension greater then one is more interesting for us, as for 
the dimension equal to one the situation is the following: if nontrivial vector 
field satisfies condition (A2), then its divergency is equal to 1 automatically. 

Further we assume that the field A has at least two times differentiable 
components. 

Let us remark that a smooth vector field with constant nonzero divergency 
cannot exist on the compact manifold without boundary. Indeed, it contradicts 
to the fact that the integral from the divergency taken over all manifold has 
to be zero. 
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Theorem 3. 2. There exists no nontrivial divergency free vector field given 
on the manifold of dimension n satisfying condition (A2). 

Denote 

Jm= Y. - V,,A*A..V,^A*-), 

ik&Km,jk¥^ik, k=l,...,m 

where the summation is taken over all subsets = ii,---fim, 1 < < n, 

from the set 1, ...,n and the set ji, is a permutation of elements of Km- 

We need the following lemma 

Lemma 3. 1. Let the field A satisfy (A2) and have a constant divergency D. 
Then the following representation holds: 

D = + const ■ Jm + -Ad(D, J 2 , Jm-i), (3.1.6) 

where A4{D, J 2 , ■■■, Jm-i) is a polynomial of k-th degree 0 < k < m, m = 

2,...,n. 

Proof of Lemma 3.1. From (A2) we get 
Ah = = A^aVigA^^Vj^A*! = 

where i^ = 1,..., n, k = 1, ..., n, m < n. Therefore 

n 

D=Y i...Vi,Wi)+ 

ii=l 

If Z1 = const, then D = D™' + const ■ Jm + Ad(D, J 2 ,..., Jm-i), the last 
polynomial is homogeneous of m-th degree with respect to VjA-^, i,j = 1,..., n. 

Example. The rather thorough (for m > 2) calculation shows that the 
following formulas are true: 

D = D^-2J2, (3.1.7) 

D = + 3 J 3 - 3 DJ 2 , (3.1.8) 

D = D^ + 4J4-4DV2 + DJ3 + 2J|. (3.1.9) 

Lemma 3. 2. For the field A satisfying (A2) the following representation 
takes place 

1 — D J 2 — J 3 J 4 — Jn = (3.1.10) 

Proof of Lemma 3.2. Consider the system, defined by the condition (A2) as 
a linear homogeneous system with respect to the variables A*. The condition 
(3.1.10) is namely the condition of existence of its nontrivial solution, that is 
the corresponding determinant of linear system is equal to zero. So the proof 


IS over. 
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Proof of Theorem 3.2. Suppose that D = 0. Then from the representation 
(3.1.6) we obtain consecutively that Jm = 0, m = 2,..., n. This fact contradicts 
to (3.1.10). 

Theorem 3. 3. Suppose that on the manifold S of dimension n there exists 
the vector field A satisfying condition (Al). If a nontrivial vector field, given 
on S and satisfying condition (A2) has a constant divergency D, then D can 
take only integer values from 1 to n. 

Proof. Equations (3.1.6) for m = 2,...,n and (3.1.10) form the system 
of n algebraic equations for n unknown variables {D, J 2 ,..., Jn)- We express 
consecutively J 2 ,..., Jn from equations (3.1.6) through the divergency D and 
substitute the result in (3.1.10). In such way we obtain the equation of n-th 
degree for D. Thus, D cannot take more then n different real values. But 
the field A, as well as its projections to the subspaces of dimension 1, ...,n — 
1, satisfy condition (A2), moreover, their divergency takes exactly n values, 
namely, 1,2, ...,n. Thus, D cannot take another values and the Theorem 3.3 
is proved. 

Remark 3.1.4. For the euclidean space A = r, where r is the radius-vector 
of point. 

Remark 3.1.5. To all appearances, the requirement of existence of the 
field satisfying (Al), is unnecessary, that is in the case where the divergency 
of the field A with the property (A2) is constant, it is equal to a natural 
number, later or equal to the dimension of space. At least the fact holds in 
the spaces of dimensions 2 and 3, as the following proposition shows. 

Proposition 3. 1. If a nontrivial vector-field given on the manifold S of 

dimension n = 2 or n = 3 and satisfying condition (A2) has the constant 
divergency D, then D can take only natural values from 1 to n. 

Proof, a) If n = 2, then from (3.1.6), (3.1.10) we have that D = — 

2 J 2 , 1 — J) -|- J 2 = 0, or — 3D + 2 = 0. The last equation has two roots 
D = l, D = 2. 

b) Analogically for n = 3 we have the system of equations D = — 

2 +2, D = D^+3Jy,—3DJ2, I—D+J 2 — JZ = 0. Therefore 11^—611-1-1114—6 = 0, 
and the roots are D = 1, D = 2, D = 3. 

3.1.2. A generalization of condition (A2). The solution to the equation (A2) 
is not unique. Suppose that we succeed in finding of finite or infinite number 
of the solutions such that Ak, k € K, K C J\f, having the additional property 

A^VjAj -|- Aj VjA(„ = /3k-^k) m / 1, /3k = const., k, m, 1 G K. (3.1.11) 

k 

Then there exists a solution of the form 

u = y~]ak(t)Ak. 
k 
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The functions ak(t) satisfy the system of equations 

= -«kW - ^kak(i) ^ ai{t), k G K. 

¥k 


(3.1.12) 


4. Two-dimensional manifold 

The section is devoted to the important for applications case of a two 
dimensional manifold. 


4.1. Transport equation on the manifold covered by one chart. Let 

be the coordinates on the manifold, and (A^,A^) be the components 
of the solution to equation (3.1.4). Thus, the equations (A2) has the form 


A^ = A 


■ + r!.(A‘)= + 2r!jA'A= + ^‘,(A^)^ 


dx^ 


dx"^ 


A2 = A^?^ + A^^ + r?i(Ai)2 + 2T\2A^A‘^ + 


dx'^ 


(4.1.1) 


(4.1.2) 


A2 
AT ’ 


then as a corollary of (4.1.1-4.1.2) we obtain the 


dx^ 

If we denote Z = 
quasilinear equation 

f)Z r)7 

+ Z-^ = Tl,{Zf + (2ri2 - Tl,){Zf + (Th - 2Tl^)Z - (4.1.3) 

which can be solved by the characteristics method: 

dZ 


dx^ 


= Tl,{zf + ( 2 rl 2 - Tl,){zf + (Til - 2 ri,)z - r?i, 


dx'^ 

dx^ 


= Z. 


If Z{x^,x‘^) is found, then from (4.1.1) we obtain the linear equation 

^ + rL(^)^)Ai = 1, (4.1.4) 

for the function A^. As A 2 = ZAi, then the solution components are found. 
If Aj = y^^A* are the physical components of the field A, then 

dXi ^/^X2^Xl . 5 iIt.i (^ 2 )^ ! - 

\/522 X\ dx^ 922 Ai 

9 X 2 , ^/^ X 2 8 X 2 \ I t^2 X A-—^2 




dx^ Xi 5x2 


+ -|- r2iA2 — 

Al 


If 2 : = ^, then the function satisfies the equation 
dz \/<7ii dz V 922 .p2 I j ^2 V 9ii 2 9'11-pi 3 


5x1 + ^^dx2 + 


+ ^212: - 


ri^z^ - = 0 , (4.1.5) 

V922 922 
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the components Ai and A 2 can be found from the equation 


<9Ai <9Ai 

dx^ dx'^ 



T\^X,z + ^Tl^z^X 

922 


1 — 




9X2 ^/^ 9X2 

9x^ 9x‘^ 



r?i^ + riiA2 = 


Here Ai = 


A 2 


(4.1.6) 


Remark 4.1.1. As we shall see, the important role plays a possibility to 
find a potential vector field satisfying (A2). In the terms of the potential 4> 
(such that A = V<I>) condition (A2) takes the form 


= (V4>,V)V$ 


or 

V^4> = V^4>Vl4> + v24>V^4>, 

V24> = V^4>Vi4> + 

If we consider the system as linear homogeneous with respect to and 
we obtain the necessary condition for the existence of the nontrivial 
potential <I> : 

Hess<k — A<I> + 1 = 0, 

where Hess<I> = V}<I>V|<k — Vf<I>V 2 <k, A<k = V}<I> + 

Remark 4.1.2. For the equation of the form 

E{x\x\A)A'^ = A^VjA\ i = l,2, 

whose particular case is (4.1.2), the equation (4.1.3) is true as well. But instead 
of (4.1.4) we obtain the equation 

DAi DAi 

^ + (Fl, + 2Fl2^ + rUzf)A^ = Eix\x\ |A|). (4.1.7) 

In the case of aerodynamic friction, where F(V) = —^iV|V|, +1 = const > 
0, we can obtain 

S(x\x^,|A|) = 1 - /ii signa(t)|A^|\/l + {Zy. 


Remark 4.1.3. If the dimension n is greater than 2, by means of intro- 
duction of a new variable Z^-i = k = 2, ...,n, acting analogously, one can 
deduce the equation (A2) to the system of n—1 equations for Zi, i = I ,..., n—1. 
In the case of euclidean metrics it will be n — 1— dimensional system of trans¬ 
port equations as well, therefore the procedure of the dimension reduction we 
can apply several times up to reduction to one equation. 
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4.1.1. Nonlinear transport equation with the Coriolis force on a two dimen¬ 
sional surface. We mean the equation 


Ut + (u, V)u = /u_L, 
or 

ui + u^V ju^ = IN.jU^, (4.1.8) 

where I = l{x^,x‘^) is a Coriolis parameter, u]_ = e^^u^, Cij is the skew- 
symmetric Levi-Civita tensor. 

We seek for a solution of the equation in the form 

u = a(t)A-|- H, (4.1.9) 

where A and H are vectors, which do non depend on time. 

In the case a'{t) = —a?‘{t) as before, the vector A satisfies condition (A2), 
the vector H is a stationary solution to (4.1.8), that is the solution to the 
equation 

i = 1,..., n. (4.1.10) 

Moreover, the vectors A and H can be compatible in some sense, that is 
have to satisfy the equation 

IN.jN = H^VjA* + NV0, i = 1,..., n. (4.1.11) 

Thus the solution of form (4.1.9) is a sum of a nonstationary solution to 
the non rotational equation (4.1.1) and a stationary solution to (4.1.8). Of 
course, it is not evident before that the condition (4.1.11) takes place. 


4.2. The case of the plane. For the euclidean metrics all Christoffel symbols 
are equal to zero, the geometrical components of the vector A are equal to 
physical ones. Therefore the equation (4.1.3) coincides with the transport 
equation (here (xi,X 2 ) = (x^,x^) are the charthesian coordinates): 


dz ^ dz 

dxi dx2 


(4.1.12) 


As well known (f.e., [17]1. the solution to equation (4.1.12) can be given implic¬ 
itly in the form z = F{x2 — zxi)., where z{ 0 , X 2 ) = F{x2) is the Cauchy datum 
on the noncharacteristic curve xi = 0. The functions Ai(xi, X 2 ), A 2 (xi, X 2 ) can 
be found from equations 

dxi ^dx2 ’ 

5A2 dA .2 

Note that the presence of singularities of a solution to (4.1.12) does not signify 
the presence of singularities of Ai,A 2 . 

The solution with linear profile of velocity of the form V = a{t)Ir B{t), 
where r is a radius-vector of point, corresponds to z = 
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It is not difficult to see that z = K = const corresponds to the velocity 
field 

Ai = Xi + (/>(X2 - Kxi), A2 = Kxi + K(I){x 2 - Kxi), (4.1.13) 

where (p is an arbitrary differentiable function. The divergency of the velocity 
field is equal to 1. 

Remark that it can be potential only if it coincides (up to the rotation with 
respect to the origin) with the linear field, depending only on one coordinate. 

Thus, we have found the method of constructing solutions with separated 
variables of the form u = a{t)A of the two-dimensional transport equation on 
the plane. 

Stress that the solution cannot be divergency free. 

Let us try to find the condition of the existence of a solution of form 

u = ai(t)Ai -I- a2(t)A2, 

where Ai is given by formula (4.1.13), and A 2 = r is a radius-vector of 
point. We can clear up, when condition (4.1.11) is true. The computation 
shows that it takes place only for 4>{^) = C^, f3i = 2 , (32 = 0, that is if 
a]; = (1 — CK)x + Cy, A^ = K{1 — CK)x + CKy. The functions ai{t), i = 1,2 
can be found from the system 

ai(t) = -a\{t) - 2ai{t)a2{t), ai{t) = -afit). 

However the obtained solution has the form u = A{t)A, A = r, where 
A'{t) = —A‘^{t), that is we have found a solution of form (4.1.2), with the 
additional condition (Al)). It is the solution with linear profile of velocity 
described before. 


4.2.1. Nonlinear transport equation on the two-dimensional sphere. For the 
sphere of radius r with the standard orthogonal coordinates x^ = (/>g] — 
TT, tt], = 0 G [0, tt] the metric tensor is the following: gn = sin^ 0, g 22 = 
r^, T^]^ = — sin0cos0, T \2 = cigO, the other Christoffel symbols are equal to 
zero. 

Thus, from equation (4.1.5) for the ratio z = ^ of physical components of 
the velocity V = (u, v) we have 

dz dO . ^ 

— = cos 0(1 -|- z ), — = sm0z. 
d(p d(p 

The equation may be solved explicitly. Namely, 


z = ± 




-1 + siffi 0T 



dr 


sinr vA + ffi sin^ T 


(4.1.14) 


with an arbitrary differential function T. Remark that this solution exists only 
in the strip | sin0| > exp(—-y), it implies K > 0. 
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Example. Consider at greater length the simplest case 'h = C = const. 
We can see from (4.1.14) that the solution is defined by the formula only 
if I sin 01 > that is in some spherical strip II surrounding the equator, 

moreover, the solution is zero on the boundary of the strip at 0* = ± arcsin 
Further, to define the meridional component v we use the equation (4.1.6), 
so that 

dv . n dv ^ 

— + sm uz— — cos 0 — = r sm Oz. 
ocp off z 

Further considerations we carry out for the nonnegative branch 2 : = \/C sin^ 0 
For the branch with the opposite sign we can act in a similar way. Therefore, 

dv . „ / „ . 9 „ —-dv cos 6 . ^ / ^ ■ 9 n —7 

7-7 + sm0v C sm 9 — 1—- , = r sm0v C sm 6 — 1. 

dO ^ycsm‘^e-l 

Solving the equation we obtain 


V = -- 


: arcsin 


\/C cos 9 \/ C sin^ 6 — 1 \J C sin^ 6 — 1 ^ ,, 

+ -Tl(0-7^(0)), 


y/c yJC - 1 sin 9 ' sin 0 

where Ti is an arbitrary differentiable function, and 


(4.1.15) 


7^ = -Intg- + 


yfC ^^tg2f-2C' + l 

, arcth _ , — . 

yfC^ 2y/Cy/C^ 


Further, for the parallel directed component we have the formula 


r \fc rns 9 1 1 

We see that the meridional component v vanishes on the boundary of the 
spherical strip II, and u remains bounded but, generally speaking, does not 
vanish. Therefore, the motion takes place inside II independently from zones 
near poles of the sphere. So the stream flows round the polar zones. 

Remark that the vector field from this example is not potential. 


5. Constructing solution to the Euler equations 

The procedure of constructing solutions to system (E1-E2) consists of the 
following steps: 


1. Find the velocity vector 


V = a(t)A(x), 


with the - smooth field A given on S satisfying (A 2 ). 


(5.1) 
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2. Supposing that we know V, find p and p from the equations (El), (E3), 
linear with respect to these functions. 


3. Choose the initial data po and po such that the compatibility condition 
Vp(t,x) = (/>(t)p(t,x)A - p(t,x)F(t,x,a(t)A,p,p) (5.2) 

holds with some function 4 ’{t). 

Remark that it may be that for a certain manifold there are no initial 
conditions po^Po and the vector field A to satisfy (5.2). 

However, as we see (section 2), there are examples, where it really holds, 
at least in the euclidean space. 

If we additionally suppose that F = F(V) and it is homogeneous with re¬ 
spect to V (the case of dry and aerodynamic friction, and of the Coriolis force 
as well), then the compatibility condition takes simpler form 

Vp(t,x) = (/>(t)p(t,x)(A - F(A)). (5.3) 

5.1. The velocity field. Integral functionals for the Euler equations. 


5.1.1. The zero right hand side. If A is found (for example by the method 
described in Section 2), then the further problem is to find the function a{t). 


Suppose that in system (E1-E3) the mass M.= pdT,, and the total en¬ 


ergy £ = Ek{t) + Ep{t) are conserved. Here the kinetic and potential energies 
are defined as Ek{t) ~ 2 J -^p(^) =-^ J pdE, respectively. 


s s 

This conservation laws always take place for the smooth solutions on the com¬ 
pact manifold. If the manifold is not compact, then it is sufficiently to require 
a rather quick vanishing of p and p at infinity with respect to the space vari¬ 
ables, to fulfill the convergence of integrals expressing the mass and energy 
and other integrals, that we mention below, as well. 

Recall that we seek for the velocity in the form (5.1). 

Introduce the functionals 


Gmit) = \j plArdS, 

s 


Q{t) = J pDdTi, D = ViA*, m > 0, 

s 

Fit) = l J p(V|A|2,V)dS. 

s 


Denote G 2 {t) := G{t). 
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5.1.2. A potential velocity field. Suppose that the field A is potential, that is 
ViA^' = VjA\ i / R A = |V|A|2. 

It is not difficult to see that here for the smooth solutions to the system 
(E1-E3) the following relations hold 


G'{t) = F{t), 


(5.1.1) 

F{t) = 2a{t)G{t), 


(5.1.2) 

G'rait) = ma{t)Gm{t), 


(5.1.3) 

F'{t) = 2Ek{t) + Q{t), 


(5.1.4) 

E'^{t) = -E'^{t) = -{^-l)a{t)Q{t), 


(5.1.5) 

Ek{t) = a‘^{t)G{t), 


(5.1.6) 

Ep{t) + a^{t)G{t) = E = const. 


(5.1.7) 

It follows that for the solutions with the velocity field of 

form 

(5.1) we 

obtain the system 



G'{t) = 2a{t)G{t), 


(5.1.8) 



(5.1.9) 

Moreover, from (5.1.3) we obtain 



Gm{t) = KmliGl{t))T, Kml = Gm{0)Gi{0)~T = cOUSt. 


(5.1.10) 

It will be more convenient to involve a new variable G{t) = 

1 

Git) 

> 0 and 

instead of (5.1.8-5.1.9) to consider the system 



a'{t) = -a^{t) + ]^Q{t)G{t), 


(5.1.11) 

G'{t) = -2a{t)G{t). 


(5.1.12) 


However this system is not closed. Below we consider the cases were there 
is a possibility to obtain it in a closed form. 

Note that the properties of the function a{t) are different from the case of 
transport equation considered in Section 3. 

Remark 5.1.1. The functional Q{t) equals to zero identically if p = 0 (so 
called pressure free gas dynamics or in the case of divergency free field A). As 
follows from (5.1.9), if a(0) < 0,then the velocity V goes to infinity, and G{t) 
goes to zero in a finite time. Taking into account the mass conservation we 
conclude that on the points where |A| = 0 the mass concentrates on the set of 
zero measure. The phenomenon really takes place in the case of the pressure 
free gas dynamics. The second possibility don’t realize as A cannot have the 
zero divergency. 
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Figure 1. 


5.1.3. A non-potential velocity field. If the field A satisfies (A2), but is not po¬ 
tential, then the properties (5.1.1-5.1.5) take place, but the expression (5.1.6) 
for the kinetic energy Ek{t) becomes more complicated. For example, in the 
two-dimensional space the following relation takes place; 

Ek{t) = a^{t)G{t) + 2 J pJh^k^{l - D) dS, (5.1.13) 

s 

with D = VjA*, J = V 2 A^ — ViA^. It is interesting that if the divergency D 
is equal to 1, then (5.1.13) coincide with (5.1.6) and we obtain once more the 
system (5.1.11 - 5.1.12). Recall that there are two possible values of constant 
divergency, namely, D = \ and D = 2 (see Theorem 3.3), the example of 
a vector field A on the plane with the divergency equal to 1 constructed in 
Subsection 2.1.1. 

5.1.4. The field A with a eonstant divergency D. Suppose once more that A 
is potential. It is easy to see that if D is constant, then 

Q{t) = ij-l)DEpit), (5.1.14) 

Ep{t) = const ■ G~' 2 ^ (t). (5.1.15) 

Thus, the system (5.1.11-5.1.12) has the closed form 

&{t) = —2a{t)G{t), a'{t) = —a^{t) -\- const ■ G 2 (5.1.16) 

The energy conservation gives 

(7-1)0 r, 

const ■ G 2 (t) + a^{t)G{t) = £. 

So for H > 0 the function G{t) cannot go to zero (a(f) to the infinity, respec¬ 
tively), therefore the mass cannot be concentrated in a point. 

The only equilibrium point is the origin. The phase portrait is presented 
on Figure 1. 
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5.2. The bounded field A with an arbitrary divergency D. If the diver¬ 
gency is not constant, then, generally speaking, we cannot obtain the closed 
system for finding a{t). However, some properties of a{t) we can get. 

It is evident that G{t) > 0. In the equilibrium points we have Q{t) = 0. It 
may exist several points of such kind on the axis a = 0. Anyway, the phase 
trajectories of (5.1.11-5.1.12) are symmetric with respect to the axis. 

Moreover, the function a{t) cannot go to the infinity. Really, from (5.1.12) 

we have that G{t) = G(0) exp ^—2 f a(T) dr^ , therefore if o(t) —oo, t t*, 

then Git'j —y oo i^Git') — y 0) t —y 

Further we use the generalization of Lemma from |19] . 


Lemma 5. 1. Let T, be a smooth riemannian manifold with the metric ten¬ 
sor gij such that \gij\ is separated from zero by a positive constant g^. If 
on the manifold there exists a function /(x) such that /(x) > 0, /(x) G 

|xp/(x) G Li(S), and 7 > 1, then 

27 n(7—1) 

f f f \ ("+2)7-" / f „ \ (n+2)7-n 

where 


"(7-1) 

_ „2((ra+2)7-n 

^i,9ij ~ y* 


27 


^re( 7 - 1 ) 


n(7-l) 
(71+2)7 —n 


+ 


27 


n(7- 1) 


27 

(71+2)7 — 71 


According to (2.3.3) the quantity S remains constant along particle trajec¬ 
tories, so that S{t,x) > sq, where sq = nun5o(x). 

We use Lemma 5.1 to obtain 

= yb/AexpSdE > ^ Jp-'dE > 

with 

1 ( M \ d’^+2)7-fi)/2 

^ = (7 _ i)2((-+2)7-n)/2 j 

and Ad is a conserved total mass. 

Therefore, if G{t) ^0, t ^ t^, then Ep(t) 00 , and we obtain the 
contradiction with the conservation of total energy. 

Remark 5.2.1 From the state equation (see Section 2.3) we have p > 
expsop^, where sq = nun5o(x). The Holder inequality gives us 


Ad^ 


< {Ll{t)y ^ j dS < exp(—So) sup H(t) J pdT, = PQEp{t), 
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with the area of the density support Q{t) and the positive constant /3o = 
exp(—So) supr 2 (t )(7 — 1). Therefore, if the support of the density is bounded, 

t 

the potential energy cannot vanish and there exists an inaccessible potential 
energy, which cannot convert to the kinetic one. Namely, 

If the support of density is not bounded, then the inaccessible energy does 
not exist, generally speaking [19] (about the treatment on the inaccessible 
potential energy see also mm)- 

5.2.1. Two-dimensional manifold. This situation is more simple. Let us in¬ 
troduce functionals 

Qm{t) = j pD^ dT,, m = 0 , 1 ,..., 
s 

where D = VjAL It is easy to see that Ep{t) = (7 — l)(5o(^)) the functional 
Q{t) introduced in Section 5.1.1 is Qi{t) in the new notation. 

Further, from (E3) for the velocity of form (5.1) we have 

Q'^{t) = a{t) I p{{A, VD^) - (7 - l)Z)-+i) dS. (5.2.1) 

From (3.1.5) and (3.1.10) we obtain (A, VD) = D — -|- 2J and J = D — 1, 

therefore (5.2.1) gives 

Q'mi't) =+ - — -)D‘^ - 3D2) dT, (5.2.2) 

J m 

s 

or 

+ + 3Qmit)-2), m E N. (5.2.3) 

For any m the coefficient at Qm+i does not vanish, therefore the system (5.1.8), 
(5.1.9), (5.2.3) is also unclosed. 

However it is interesting that 

sign( 3 („(t) = -signa(t) 

for m = 2fc -|- 1, /c = 0,1,... , m < 8(7 — 1). It follows from (5.2.2), as the 
expression 

+ -3D+ 2 

m 

is positive and separated from zero for m < 8(7 — 1 ). 

In particular. 


Qi{t) = -a{t)N{t) 


( 5 . 2 . 4 ) 
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Figure 2. 


where N{t) 


J p{'yD‘^ — 3D + 2 )pdT,, moreover 
s 


N{t) > aoEp{t), ao 


(87-9)(7-1) . 9 

- ^ -' ^^ 8 - 


From (5.1.8) and (5.2.4) we have —^7 = —77 < 0. Therefore there exists a 

dCj 

decreasing function Q\ = Q\{G) and the system (5.1.8), (5.1.9) can be written 
as 


G'{t) = 2a{t)G{t), a'{t) = + - 


IQi(G) 


G 


(5.2.5) 


Let us study the equilibriums of (5.2.5). Remark that at the phase plane 
(G, a) the trajectories are symmetric with respect to the axis a = 0. The 
second coordinate of the equilibriums is a = 0, the first one is G satisfying 
the equation Qi{G) = 0. Recall that if the root of the equation exists, it 
is unique. If the root does not exist, the first coordinate of the equilibrium 
goes to the infinity, in the case if the field A is unbounded. Otherwise, if 
G{t) < G+ = const, it signifies that there exists no globally in time smooth 
solution with the properties we try to construct. If we use once more the 
variable G = G~^, then we get the system 


G’{t) = -2a{t)G{t), a'{t) = -a\t) + ^Qi{G)G, 


and the equilibrium point goes from the infinity to the origin. The phase por¬ 
trait is similar to the one from Fig.l. If the field A is bounded, the restriction 
for G has the form G > The situation is presented on Fig.2. 

If the manifold is compact, then G{t) < —A4max|Ap (if we suppose A to 

be smooth), therefore there is no equilibrium in the infinity. Note that Qi{t) = 
0 for the constant solution p = po, Y = 0 , p = pq. Here the coordinates of 
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Figure 3. 



Figure 4. 


the only equilibrium are (Go, 0 ), with Go 



Ap dS = const > 0 , it is a 


center (see Fig.3). 

If the manifold is not compact, it may exist two equilibriums - in the finite 
point (Go,0), where Go is a root of the equation Qi{G) = 0 (a center) and in 
the infinite point (G = +oo, a = 0), (the point (G = 0, a = 0) on the plain 
(G, a) respectively). The situation is presented on Fig.4. The infinite point 
arises only if A is unbounded. 

Remark 5.2.2 If we suppose that there exists the limit function Qooit) = 


lim 

m^oo 


dTi, then from (5.2.3) we obtain that the equation 


s 


Qoo{t) = -(7 - l)a{t)Qoo{t) 
holds. Together with (5.1.8) it gives that 

'y — 1 

Qoo{t) = const ■ G ~{t). 
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Let US stress that it is not clear whether the limit exists in the nontrivial case 
Qoo / 0. 

5.2.2. Appearance of singularities of smooth solutions. Suppose, for example, 
that the smooth field A is bounded together with its divergency D. Generally 
speaking, the divergency changes the sign over the manifold, that is inf D < 0. 
From (5.1.4) we obtain 



(5.2.6) 

where A^ = sup Ap, T)__ = inf D\. 


Moreover, 


F{t) < 2Ek{t)G{t) < MEK\. 

(5.2.7) 

It follows from (5.2.6) that if we choose 


F( 0 ) > A+V( 7 - 1 )H_MF, 

(5.2.8) 


then F{t) will grow unboundedly, it results in a contradiction with (5.2.7). 

As the mass does not concentrate on sets of zero measure, that F{t) will 
grow owing to the gradient of density. Really, in the case of a potential field 

A we have F{t) = /(pV, A)dS = -/V(pV)^(iS = -/(Vp, A)^dS - 
s s s 

f The last integral is bounded by the constant, therefore, the grow 

s 

of F(t) corresponds to the growth of |(V/j, A)| with respect to time. 

Remark that one can ask whether the condition (5.2.8) is possible? The 
comparison of (5.2.7) and (5.2.8) gives the necessary condition for this possi¬ 
bility, namely, 

7-1 

5.3. The Euler equations with the right hand side of special form. 

Let A satisfy (A2) and the velocity field be of form (5.1). Suppose that for 
F(V) the following condition holds 

K 

(F(V), A) = ^Afc(^,a(^))|Ar^ (5.3.1) 

k=l 

k G M, Sk > 0. Here Xk{t,a{t)) are known functions of a{t) and t. 

Remark 5.3.1 The condition (5.3.1) takes place, for example, in the case 
of dry and aerodynamic friction. 

Recall that for the dry friction 

F(V) = -/i(t,x)V. 

Suppose fi{t,x) = ii{t). Then A = 1, si = 2, Xi{t,a{t)) = —iJ,{t)a{t). 
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For the aerodynamic friction 

F(V) = -/xi(t,x)V|V|. 

If then K = 1, si = 3, Xi{t,a{t)) = —/ii(t)a(t)|a(t)|. 

Remark 5.3.2 Also (5.3.1) holds for the Navier-Stockes equation in the 


Euclidean space, where E* = V-^Tj, with the tensor 


r; = a{p){ViViV^ - VjV^) + (5.3.2) 

if A coincides with the radius-vector r ((F(V) = 0 in the situation). 

Thus, from (5.1.1), (5.1.2) we obtain similar to (5.1.3) that 

K 

F'{t) = {2a{t)G{t)y = 2Ek{t) + Q{t) + ^ Afc(t, a{t))Gs, {t), 


k=l 


and 




k=l 


_a. 

According to (5.1.10) we have Gs^{t) = Kg^G^ (t) = Ks^,G~ (t), with some 
constants depending on the initial data. Thus, 

a'(t) = -a^(t) + a{t))G^-\t) + 


k=l 


Introduce for the convenience (similar to Section 5.1.1) the new variable 
G{t) = G~^{t), and obtain the system for (a(t), G{t), Q), n G J\f: 

. K ,, 1 - 

a'{t) = -a^{t) + Ks^ 2 Xk{t, a{t))G^~^ {t) + -Q{t)G{t), (5.3.3) 

^ k=l 

G'{t) = -2a{t)G{t). (5.3.4) 

Generally speaking, it is unclosed, but in the case of constant divergency 
VjA* = D, as in the case of zero right hand side it is possible to express Q{t) 
through G{t) and obtain the closed system for {a{t),G(t)). Namely, 

1 K 

, o i X— > ~ ~ ( 7 — 1 )£) + 2 

a {t) = -a^{t) + 2 X] Ks^^ 2 Xk{t, a{t))G^~^ (t) + KG 2 (t), (5.3.5) 

^ k=l 

G'{t) = -2a{t)G{t). 

For the dry friction with constant friction coefficients p, the system (5.3.5 
- 5.3.4) takes the form 


. n 1 ~ ~ (7-1)D + 2 

a {t) = -a^{t) - -pKs^ 2 a{t) -f KG 2 (t), 

G'{t) = -2a{t)G{t). 


(5.3.6) 



42 


OLGA S. ROZANOVA 



Figure 5. 



Figure 6. 

For the aerodynamic friction with constant coefficients we have 

a{t) = -a^{t)--^ilKs^2a{t)\a{t)\G ■^{t) + KG 2 {t), (5.3.7) 

G'{t) = -2a{t)G{t). 

The phase portraits of systems (5.2.6 - 5.2.4) and (5.2.7 - 5.2.4) are pre¬ 
sented on Figures 5 and 6, respectively. 

In the case of dry friction there are two equilibriums: at the origin (stable) 
and at the point (0, (unstable). 

In the case of aerodynamic friction (and for any friction term with a > 0 
as well (see Section 1.1)) the only equilibrium is at the origin, moreover, it is 
stable. 


6. Further discussion 

There is a lot of interesting problems connecting with this work that would 
be a subject of further investigations. Let us mention some of them. 
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1. Euler equations on the two-dimensional sphere 

a) to construct, if it is possible, a potential smooth vector field satisfying 
(A2). 

b) to find a possibility to close the system (5.1.11 - 5.1.12) and to investigate 
the behaviour of the function a{t). 

2. Interior solutions 

It is rather easy to show that if in the Euclidean space the vector field 
A satisfying (A2) differs from the radius-vector up to the function from the 
Sobolev space, then the corresponding solution is interior. 

Which of solutions with velocity field (5.1) are interior? At least is it true 
that any solution where A has a constant divergency is interior? 

3. Nonviscous solutions of systems with viscous term 

It is easy to see that in the Euclidean space the solutions with linear profile 
of velocity ”do not feel” the viscous term (5.3.2) in the Navier-Stokes system. 
Therefore the solutions solve Navier-Stokes as well. In particular, for the 
viscous, but pressure free case we can construct the solution with the mass 
concentrated in the point. 

It is interesting to find the solution of such kind for any manifolds. What is 
a condition for the metric tensor for the possibility of existence of ’’nonviscous 
solutions” ? 

Further, if a solution with linear profile of velocity is interior for the Euler 
system, can we guarantee that it will be interior for the corresponding Navier- 
Stokes system? 

Acknowledgments. The work is partially supported by the Russian Foun¬ 
dation for Basic Researches Award No. 03-02-16263. 

References 

[1] T. Makino, S.Ukai, S.Kawashima (1986): Sur la solutions a support compact de 
I’equations d’Euler compressible. Japan J. Appl. Math., 3, 249-257. 

[2] B.A.Dubrovin, S.P.Novikov, A.T.Fomenko (1992): Modem geometry - methods and ap¬ 
plications. Part I. The geometry of surfaces, transformation groups, and fields. Gradnate 
Texts in Mathematics. 93. New York etc.: Springer-Verlag. 

[3] T.Kato (1975): The Canchy problem for quasilinear symmetric hyperbolic systems. 
Arch. Ration. Math. Anal.,58, 181-205. 

[4] A.I.Volpert, S.I.Khudiaev (1972): On a Cauchy problem for composite systems of non¬ 
linear equations. Mat.Sbornik 87, N4, 504-528. 

[5] L.Carding (1963)'.Probleme de Cauchy pour les systemes quasilineaire d’ordre un stricte- 
ment hyperboliques. in Les equations anx derivees partielles, Coloqnes Internatiananx 
de CNRS, vol.117, 33-40, Paris. 

[6] P.D.Lax.(1973) Pbyperbolic systems of conservation laws and the mathematical theory of 
shock waves. SIAM Reg.Conf.Lecture, Nil, Philadelphia. 

[7] A.Majda(1984): Compressible fluid flow and systems of conservation laws in several 
space variables. Applied Math.Sci.V.53. Berlin: Springer. 

[8] M.Crandall, P.Souganidis (1986):Nonlinear Analysis 10, N5. 



44 


OLGA S. ROZANOVA 


[9] J.-Y. Chemin (1990): Dynamique des gaz d masse totale finie Asymptotic Analysis. 3, 
215-220. 

[10] O.S.Rozanova (1998): Generation of singularities of compactly supported solutions of 
the Euler equations on a rotated plane. Differ. Equat., 34(8), 1118-1123. 

[11] Z. Xin(1998): Blowup of smooth solutions to the compressible Navier-Stokes equation 
with compact density. Comm.Pure Appl.Math., 51, 229-240 

[12] O.S.Rozanova (1998): Energy balance in a model of the dynamics of a two-dimensional 
baroclinic atmosphere. Izvestiya, Atmospheric and Oceanic Physics, 34(2), 186-196. 

[13] V.A.Gordin (2000): Mathematical Problems and Methods of hydrodynamic Weather 
Eorecasting, London: Taylor & Francis Group. 

[14] O.S.Rozanova (2002): On classes of globally smooth solutions to the Euler equations in 
several dimensions. LANL e-print math.AP/0203230 

[15] M.Grassin (1998): Global smooth solutions to Euler equations for a perfect gas. Indiana 
University Math.J.,47(4), 1397-1432. 

[16] D.Serre (1997): Solutions classiques globales des equations d’Euler pour un fiuide parfait 
compressible. Annales de I’lnstitut Fourier, 47, 139-153. 

[17] G.B.Whitham (1999): Linear and nonlinear waves. Pure and Applied Mathematics. A 
Wiley-Interscience Series of Texts, Monographs, and Tracts. New York, NY: Wiley. 

[18] R.Bellman (1960): Introduction to matrix analysis. Mograw-Hill Book company, inc. 
New York, etc. 

[19] J.-Y. Chemin (1990) Dynamique des gaz d masse totale finie Asymptotic Analysis, 3, 
215-220. 

[20] O.S.Rozanova (1999) Blow-up of solutions in system of atmosphere dynamics, In: Hy¬ 
perbolic Problems: Theory, Numerics, Applications, 

ISNM(Int.Ser.Numer.Math.)130, Vol. II, pp.793-801, Birkhauser/Basel, Switzerland. 

Mechanics and Mathematics Faculty, Moscow State University, Leninskie 

Gory 1, Moscow, Russia. 

E-mail address'. rozanova@mech.math.msu.su 


